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Lattice Gross-Neveu model with domain-wall fermions* 
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We investigate the two-dimensional lattice Gross-Neveu model, using the domain-wall fermion formulation, as 
a toy model of lattice QCD. We study features of the phase diagram related to the mechanism of chiral symmetry 
restoration, and find that the parity-broken phase (Aoki phase) exists for finite extent in the extra dimension 
{Ns). We also find that 0{a) scaling violation terms vanishes in the limit of Ns oo. 



1. Introduction 

Defining chiral fermions on the lattice has been 
one of long-standing problems in lattice field theo- 
ries. Several years ago, domain-wall fermion [|l],|j 
has been proposed as a new formulation of lattice 
chiral fermion. This formulation considers Wilson 
fermion in D -f 1 dimensions with the free bound- 
ary condition in the extra dimension of a size , 
or equivalently, TV^-flavored Wilson fermions with 
flavor mixing. In the limit Ns oo, the spec- 
trum of free domain- wall fermion (DWF) contains 
massless modes at the edges in the extra dimen- 
sion. The massless modes are stable under per- 
turbation from weak gauge fields. 

While there have been numerical simulations 
of lattice QCD with DWF (DWQCD), some non- 
perturbative issues, in particular existence of the 
parity-broken phase (Aoki phase) (|] and necessity 
of fine tuning of couplings to restore chiral sym- 
metry for finite Ns, have not been clarified. In 
order to answer these questions, we investigate 
the two-dimensional lattice Gross-Neveu model 
using the DWF formalism (DWGN) g as a test 
of DWQCD. 

2. Action and Effective potential 

We propose the following action for DWGN: 

S = Sfree + X! + *75n(n)} q{n 



+a 



.(1) 



The auxiliary fields are related to fermion con- 

2 

densates as <T{n) — irif — j^q{n)q{n) and n(n) = 

2 

— ^q{n)i^c^q{n). The important point in this ac- 
tion is that the interaction terms are constructed 
from the edge state: q{n) = Puipin, s = 1) + 
PLij[n,s = Ns). 

We attempt to solve DWGN analytically in the 
large N limit. The effective potential can be cal- 
culated by the technique of the propagator ma- 
trixl^,^, which yields 
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/ = /in [Fa^{(7^ + tf) + Gaa + H] 
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(2) 
(3) 
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where = Ili°/a (t-Tr)'^ ■ '^^'^ explicit form of the 
functions F, G and H are described in Ref. Q . Let 
us note that the term Gacr explicitly breaks chiral 
symmetry. 

3. Existence of Aoki phase 

The phase diagram of the Wilson fermion ac- 
tion has a region of spontaneously broken parity- 
flavor symmetry (Aoki phase), which plays an 
important role for controlling restoration of chi- 
ral symmetry in the continuum limit Since 
DWF is an extension of the Wilson fermion for- 
malism, we wish to examine if Aoki phase exists 




for DWGN. In this section, we set the coupling 
constants as g'^ = = g"^ . 

Let us first consider the case oi Ns ~ oo. Since 
F and H dominates over the term G in this Umit, 
the effective potential becomes 

/(a,n)= I \n[Fa^{a'' + I{^) + H] . (4) 

The 0(a) term, which breaks chiral symmetry, 
is absent, and hence the model has exact chiral 
symmetry. Thus pion becomes massless without 
fine tuning even for finite lattice spacings and for 
arbitrary strong coupling. 

Next we consider the case of finite Ns- It is 
expected that Aoki phase exists in this case since 
the Ns = I DWGN is equivalent to GN model 
with Wilson fermion0|. 

The solution of the gap equation marking the 
phase boundary is illustrated for several values of 
M with a fixed size TV^ = 2 in Fig.|]. As is sum- 
marized in a schematic diagram in Fig.||, we find 
that (i) Aoki phase exists inside the boundary 
and this boundary forms cusps toward 5^ — > 0, 
(ii) the Aoki phase always appears in the m/ > 
region for even Ns (in the conventional choice of 
sign in domain- wall literatures, this corresponds 
to TO/ < 0), (iii) pion becomes massless on the 
boundary. The last point means that pion mass 



does not vanish at m/ = for finite lattice spac- 
ing. Hence a fine tuning is needed to obtain mass- 
less pion. 

In Fig.|| we show the Ns dependence of the 
phase boundary for Ng =2,4 and 6 at M = 0.9. 
We find that (i) the Aoki phase shrinks exponen- 
tially with Ns, (ii) the first "finger" on the left 
approaches mf — Q, while other "fingers" move 
toward mf = 00, (iii) the critical value gc (see 
Fig.^) , below which the cusp structure ( "fingers" ) 
appears, increases exponentially with Ns, in con- 
trast to the case of Wilson fermion^. These fea- 
tures mean that the area of the normal phase be- 
comes wide with increasing Ns. 

As seen from above, DWF for finite Ns repre- 
sents an improved Wilson fermion. 

4. 0(a) scaling violation (a 0) 

We study the mechanism of chiral symmetry 
restoration in the continuum limit. In particular, 
we wish to understand scaling violation for finite 
and infinite Ns- 

The effective potential in the continuum limit 
is obtained after some calculation: 

K,/ = -man + ^ (4 + ^1) In • (5) 
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Figure 3. N s dependence of the phase boundary 
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Figure 4. Scahng of a/ A as a function of aA 



(1 - M)^= 

Here gb. = hi{<^ - — } and Ili^ = /a/II 

with f A f = — M{2 — which is the normal- 

^ 1 - (1 - Mf^^^ ' 
ization factor of the edge state "g(n)" for finite 
Ns- In order to obtain (H), which agrees with the 
continuum result, we need to impose the foUowing 
scahng relations: 
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These relations are the same with those found 
for the Wilson fermion case 0. Therefore a fine 
tuning is needed to restore chiral symmetry. 

Figure ^ shows cr/A as a function of a A for 
M — 0.9, using the Wilson-like scaling relations 
in (^^). We find that 0{a) scaling violation is 
large at Ns = 2,3. However, the magnitude of 
0{a) scaling violation diminishes exponentially as 
Ng increases. In fact the scaling curve almost 
exactly follows the 0{a?) behavior for A^s > 8. 

5. Conclusions 

We have investigated the two-dimensional 
DWGN model in detail as a test of DWQCD. 



When the size of the extra dimension Ng = oo, 
the model has chiral symmetry for finite lattice 
spacing, and Aoki phase does not exist. On the 
other hand, in the case of finite Ns, Aoki phase 
does exist and a fine tuning is needed to restore 
chiral symmetry in the continuum limit. How- 
ever, the 0{a) scaling violation that gives rise to 
this behavior vanishes exponentially fast as Ng is 
increased so that it is negligible in practice for 
Ns = 0(10). 

While the GN model does not have gauge fields 
and quantum fluctuations are absent in the large 
N limit, it is expected that the results obtained 
in this work provide instructive and systematic 
information for DWQCD simulations. 

The authors are JSPS Research Fellows. 
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